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Abstract 

We study the production of cc pairs in nuclear reactions at SPS energies 
within the covariant transport approach HSD. The production of cc is treated 
perturbatively employing experimental cross sections while the interactions of cc 
pairs with baryons are included by conventional cascade-type two-body collisions. 
Adopting 6 mb for the cc-baryon cross sections the data on J /'^ suppression in 
p + A reactions are reproduced in line with calculations based on the Glauber 
model. Additionally the dissociation of the cc pairs by strings is included in a 
purely geometrical way. We find good agreement with experimental data from 
the NA38 and NA50 collaboration with an estimate for the string radius of 
Rs PS 0.2- 0.3 /m. 
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1 Introduction 



The J/\l/ suppression has been proposed as one of the signals for the quark-gluon plasma 
(QGP) which is expected to be formed in nucleus-nucleus collisions at sufficiently high 
energies |jl|]. This suggestion has stimulated a number of experiments, which indeed 
have observed a significant reduction of the scaled J /'^ yield when going from proton- 
nucleus to nucleus- nucleus reactions [0]. Especially the NA50 experiment has reported 
an abrupt decrease in J/^> production in Pb + Pb collisions at 158 GeV per nucleon 
^, 1^ in going from peripheral to central collisions. 

A lot of theoretical effort has been spent to understand the experimental results 
(see for a recent review). Beside the suggestion of a possible formation of a QGP 
0, H, 3 various scenarios based on J/\Ef absorption by hadrons have also been proposed 
, [TT], |T^, |T3|, |I^, |TB[ . A part of the suppression can be explained by J / ^/^-absorption on 



the surrounding nucleons |T0[, additional absorption might be attributed to 'comovers' 



'mesons') being produced as secondaries 0, O. It has been shown in within the 



microscopic covariant transport approach HSD (Hadron String Dynamics) that 
the observed suppression of the J/\& yield in nuclear collisions is consistent with such 
an hadronic absorption scenario. However, the 'comover' models are still a matter 
of debate: The employed absorption cross section of cc pairs on comovers of 3 mb is 
treated as a free parameter in order to explain the data. On the other hand, this cross 
section might be overestimated considerably . 



In our present work we therefore consider an alternative mechanism for J/\I' pro- 
duction in heavy ion collisions. We will focus on the effect of cc dissociation in the 
prehadronic phase and not on absorption by comovers. This is motivated by the fact, 
that the very early collision phase is not described by hadrons but by highly excited 
strings, which in the HSD model are created by momentum transfer among target and 
projectile nucleons . Their production and decay is treated within the FRITIOF model 
|T8| and describes the first few fm/c of the collision. The cc state, most likely as a color 
octet, is also produced at the earliest state of the reaction in a hard collision among the 
nucleons via gluon-gluon fusion. Hence it is natural to ask what happens if this colored 
state moves into such a temporary environment of strings. As the string carries a lot of 
internal energy (to produce the later secondaries) in a small and localized space-time 
volume the quarkonia state might get completely dissociated by the intense color elec- 
tric field inside a single string. In this respect Loh et. al. have investigated the 
possible J/\E' dissociation in a color electric fiux tube in a semiclassical model based 
on the Friedberg-Lee color dielectric Lagrangian. They find that the cc state becomes 
dissociated rather immediately on a timescale less than 1 fm/c. In a heavy ion collision, 
especially for the more heavy systems, the effective region (or volume) of all the strings 
being produced is expected to be large so that the strings might become rather closely 
packed. If the cc states will get dissociated by the numerous and individual strings this 
will lead to an additional suppression of J/'^ in the early phase prior to hadronization. 
To explore this intuitive idea we include the effect of cc dissociation by strings into the 
HSD model and study the overall J/'^ production and dissociation dynamically. 

The HSD model is briefly presented in section 2, while in section 3 we describe the 
cc production process and show details of the string evolution in the HSD model. We 
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then address the dissociation process of the cc states on baryons and the chromoelectric 
flux tube of the strings (section 4). We conclude with a summary of our investigations. 



2 The covariant transport approach 



In this work we perform our analysis along the line of the HSD approach |T^ in the 
cascade modus which is based on a coupled set of covariant transport equations for the 
phase-space distributions fh{x,p) of hadron h W^, i.e. 



J2 [ d2d3d4 . . . lGtG],2^34. + p5 - - PJ ■ ■ ■) 

2,3,4... •' 

- fl{x,p)f2{x,P2)f3{x,P3)h{x,P4)} ■■■ ■ (1) 

Here [C^ G]i2^3i..S'^{pi + P2 — P3 — Pa • • •) is the 'transition rate' for the process 1 + 2 —> 
3 + 4 + . . ., while the phase-space factors 

f\{x,p) = l±Mx,p) (2) 

are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depend- 
ing on the type of hadron in the final/initial channel. The dots in eq. (|ID stand 
for further contributions in the collision term with more than two hadrons in the fi- 
nal/initial channels. The transport approach (|I]) is fully specified by the transition rates 
G'^G6'^{. . .) in the collision term, that describes the scattering and hadron production 
and absorption rates. This transport approach was found to describe reasonably well 
hadronic as well as dilepton data from SIS to SPS energies ||16|, ^ . 



In the present approach we propagate explicitly - apart from the baryons (cf. [|I6|) - 
pions, kaons, r^'s, //"s, the 1~ vector mesons p, u, (j) and K*^s as well as the axial vector 
meson ai. The high energy hadron-hadron collisions are described by the FRITIOF 



model |T8| resulting in two excited strings. The dynamical evolution of the strings is 
now included explicitly. A string is characterized by the leading quark and diquark (or 
antiquark in the case of a mesonic string) and by the energy stored in between. The 
fragmentation of the strings into hadrons starts after the formation time, which is set 
to Tp = 0.8 fm/c (see Fig. |lD. Tp controls the baryon and meson rapidity distribution 
dN/dy in comparison to experimental data. The length of the strings is given by the 
(center of mass) collision time to of the hadrons, the formation time and the velocity 
of the leading quarks/di quarks /3i = « = 1, 2 

\I32- l3l\-{t-to) , {t-to)<Tp 

l{t) = { (^|/?2-/?l|-2-yi- -(t-to) , TP <{t-U)<tmax (3) 

, (t — to) > tmax-, 
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where tmax = + 2Tj7/(y4— \f32 — is the time when the string is completely 
hadronized. The radius of the flux tube is an unknown parameter, which will be 
of particular significance for the cc dissociation by strings, which we will address in 
section 4. The cross section of the high energy secondary interactions of the leading 
quarks/diquarks are treated within a simple additive quark model a{q—B) = 1/3 a{B— 
B) and a{qq - B) = 2/3 a{B - B). 

In our simulation the production and decay of strings dominates the very early 
collision phase before the subsequent hadronic state of matter is formed. In Fig. |^ 
a characteristic representation of the hadrons and the strings during the high density 
phase in a central Pb + Pb collision at 160 AGeV is shown. As pointed out already in 
the introduction several hundred strings are formed during a central Pb-Pb collision at 
SPS energies. It turns out that most of them are in fact rather short due to secondary 
interactions of the leading quarks/diquarks. In Fig. ^ the number of strings and the 
averaged string length in a central Pb-Pb collision at IGOGeV is shown as a function 
of the center of mass time. 

After the hadronization many low energy rescattering processes are included. The 



meson-baryon and baryon-baryon collisions are treated similar to the BUU model [21 



As meson-meson channels we include the reactions vrvr — p, vrvr — KK^ vrp — 0, vrp ^ 
ai as well as the time reversed reactions using Breit-Wigner cross sections with param- 
eters from the literature |2^ and exploiting detailed balance. For the present analysis 
the low energy rescattering processes do not play an essential role. Within the phi- 
losophy outlined in the introduction the production and absorption of the cc states 
happens in the first few fm/c of the collision before secondary particles are produced. 



3 Charmonium production 

Since the probability of producing initially a cc pair is very small, a perturbative 



approach is used for technical reasons as described in |]I4[. Whenever two nucleons 
collide a cc pair is produced with a probability factor W , which is given by the ratio 
of the J/\E' to NN cross section at a center-of-mass energy a/s of the baryon-baryon 
collision. 

The parameterization used for the J/\I' cross section, the rapidity and pt distribution 
of the cc states is taken as in , which fits reasonable well the experimental data. 

It is known from experiment that the Drell-Yan cross section as a hard process 
scales with [Ap x At) 0, |^, This is in contrast to most of the soft hadronic 

observables, like e.g. the pion multiplicity, which scale like Ap + Ap. This indicates a 
strong difference in the production of hard and soft processes in heavy ion collisions. 
Since the production scheme for cc is similar to the production of Drell-Yan pairs the 
total initial cc production cross section also scales with Ap x Ap. It is important 
to stress that only with this stringent assumption the ratio of the J/\E' to the Drell- 
Yan cross section is a direct measure for the J/\Ef suppression. We want to note that 
the Ap X Ap behavior of the initial cc production cross section is still a matter of 
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debate. As pointed out by Frankel and Frati and also by Tai An et al. the 
strong J suppression in central Pb+Pb collisions might be explained by a nonlinear 
scaling of the initial cc production cross section due to initial state interactions. Such 
initial state interactions should also diminish the Drell-Yan cross section, which, on 
the other side, is not seen experimentally (for a comprehensive discussion on this issue 
we refer to [^). In any case, in our approach we stay with the same scaling behavior 
and investigate the subsequent absorption of cc pairs on baryons and strings. For that 
reason we emphasize that the ApX At scaling of the initial cc production is an input of 
our calculation and not a result. To implement this scaling we separate the production 
of the hard and soft processes: The space-time production vertices of the cc pairs are 
calculated before every run neglecting the soft processes. After that we follow the 
motion of the Charmonium pairs within the full hadronic background by propagating 
it as a free particle. Again, only with this concept the Ap x A^ dependence of the cc 
production can be reproduced in a hadronic transport simulation. 



4 cc dissociation 

The cc state as a rather heavy hadronic particle and a result of a hard process is formed 
immediately in comparison to the soft particle production from string fragmentation. 
Thus first of all the cc states move not in a hadronic environment but in an environment 
of color electric strings of 'wounded' nucleons (see Fig. |^). As motivated in the intro- 



duction and in ref . |]T9|] we now assume that a cc state immediately dissociates whenever 
it moves into the region of the color electric field of a string. In this sense strings are 
completely black for cc states. It was found in []T^ that the additional force acting on 
the charm quarks is given by 2 x a ~ 2 GeV/ fm, where a denotes the phenomenolog- 
ical string constant of a chromoelectric flux, which is sufficient to immediately break 
up a cc state. One can also argue that the field energy density contained in a string 
is given by a/{nR'g). For Rs ~ 0.3 /m one accordingly has a local high color electric 
energy density of ~ AGeV/ fm^, which substantially screens the binding potential of 
the Charmonium state 

For practical reasons the dissociation by a string is modeled when the center of mass 
of the cc state is located inside the string. The length of the strings is given by eq. (0), 
while the string radius Rs is an unknown parameter. Absorption by strings spanned 
between the parent particles of the cc pair are excluded, since this effect already is 
included in the production cross section. 

Additionally, the cc pair may be destroyed by collisions with incoming baryons. For 
the collisions with baryons we use the minimum distance concept described in Ref. . 
For the actual cross sections employed we assume that the cc pair initially is produced 
in a color-octet state and immediately picks up a soft gluon to form a color neutral 
cc — g Fock state 0] (color dipole). This extended configuration in space is assumed 
to have a 6 mb absorption cross section in collisions with baryons (cc + B Ac + D) 
as in Refs. |10|, 0] during the lifetime r of the cc — g state, for which we adopt 
r = 0.3 fm/c as suggested by Kharzeev [^. One also has to specify the absorption 
cross sections of the formed resonances J/\E' on baryons. For simplicity we use 3 mb 
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following Refs. Wl p 



5 Results 

In proton-nucleus collision the cc dissociation by strings is practically negligible. The 
absorption is dominated by cc-baryon collisions so that the proton nucleus data allow 
to fix this absorption cross section. In Fig. ^we show our results for the J/^ survival 
probability S*"^/* using 6 mb for the absorption cross section of the cc-pairs on nucleons 
in comparison to the data The experimental 'survival probabilities' 5"//* in this 
figure as well as in the following comparisons are defined by the ratio of experimental 
J/^ to Drell-Yan cross sections as 

yCTAB 1 2.9-4.5 GeVy / ^ CFpd J 

where A and B denote the target and projectile mass while a'H^ and stand for 
the J /'^ and Drell-Yan cross sections from AB collisions, respectively, and B^^ is the 
branching ratio of J/\E' to dimuons. The theoretical ratio is defined as 

where iVj/^ is the multiplicity of initially produced J/\E''s while Mj/q, is the multiplicity 
of J/\&'s that survive the hadronic final state interactions. 

In Fig. ^ results are shown for two different string radii Rs = 0.1 /m and Rs = 
0.4/m. The difference of the curves is rather small, which indicates the small effect of 
dissociation by strings. For p + U and Rs = 0.4/m only 2% of the J/\E' are absorbed 
by strings. This is completely different for heavy ion collisions, where the absorption 
on strings becomes a much more important effect. 

To compare our results for S + U and Pb + Pb to the NA38 and NA50 data, the 
experimental trigger conditions must be included. In these experiments only events 
are recorded with a fi~^fi~ pair of invariant mass M > 1.5 GeV. This trigger condition 
is obtained by 



/n /V 
^^^^w E^^rw, (7) 



where Wf^'^{b) are the weights for produced /i/i pairs within the experimental cuts. A^o 
is a normalization factor to adjust to the experimental number of events. We take the 
same weight factors W-^'^{b) as in which reasonably reproduce the experimental Et 



distribution as shown in fl^ . 

Qualitatively the results are not changed by this procedure. In Fig. ^ our results 
are shown for S + U and Pb + Pb as a function of the transverse energy and for four 
different string radii. A strong dependence on the string radius Rs is observed and 
Rs = 0.2 fm gives the best fit to experimental data 0, 0. With this string radius 40% 
of the absorbed J/^^'s are dissociated by strings in central collisions of Pb -|- Pb as 
shown in Tab. |l|. 
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Rs[fm] 


Atot 


Ab 


As 


0.1 


0.65 


0.53 


0.12 


0.2 


0.74 


0.44 


0.30 


0.3 


0.82 


0.36 


0.46 


0.4 


0.88 


0.29 


0.59 



Table 1: Total J absorption probability Atot = 1 — S^l^^^, absorption by baryons Ab 
and by strings As for different string radii in a central Pb + Pb collision as obtained 
in the HSD approach. 



6 Summary 

In this work we have addressed the question of cc absorption in relativistic heavy ion 
collisions in the microscopic transport approach HSD. Of particular interest was the 
effect of the prehadronic phase, which in the HSD model is described by independent 
strings. In principle the comover absorption scenario could be included in addition. 
However, since most of the cc pairs are already absorbed by strings and baryons in 
the early stage of the reaction, the effect of comovers is expected to be much less than 
found in ref . M . 



The dissociation mechanism by strings was treated in a simple geometrical picture. 
Whenever a cc pair moves into a string it dissociates and forms DD mesons. Strings 
are completely 'black' for cc states. This stringent condition is a first ansatz to inves- 
tigate the effect of the cc dissociation by strings. Better theoretical foundations of the 
dissociation mechanism of charmonium {J/'^ and ip') and preresonant states {ccg) are 
definitely needed to improve the understanding of this geometrical picture. 

We found that the absorption by strings plays is an important effect in the first 
few fm/c of the collision phase before secondary particles are produced. Adopting a 
string radius of 0.2/m we got a qualitative agreement with p+A data and the NA38 
and NA50 data. This radius seems to be rather small, but for two reasons it should 
be seen as a lower bound for Rs- First of all, we have assumed that the strings are 
completely black for cc states, and secondly it dissociates whenever it moves into the 
region of the color electric field of a string. With the requirement that the total cc 
state should be inside the string when the dissociation process starts, one should add 
the cc radius to our value of Rs- This would give a string radius of Rs ~ 0.4 — 0.5 fm. 

With our results we can estimate the fraction of the reaction volume which is filled 
by the flux tubes of the strings. In the very central region with 1 fm thickness along 
the longitudinal direction in the cms frame, we get Ns ~ 420 strings in a central 
Pb + Pb collision generated within the HSD approach. The averaged string length 
< / >~ 0.5/m is rather small due to secondary interaction of the leading quarks. Only 
a small amount of strings reaches a length above 1.2 fm. With these results we get a 
string volume of 



Vs = nRl-<l>-Ns^<. 3 



65/m3, Rs = 0.3/m 
115/m3, Rs = 0.4/m 
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The total volume in this central cylinder is given by 



Vt = ttR'p, ■ Ifm ^ 130/m^ (9) 

where Rpi, = 6.6 fm is the radius of the Pb nuclei. Comparing Vg and Vt illustrates 
that 50-75% of the central volume in an ultrarelativistic heavy ion collision is filled by 
strings in the first few fm/c assuming a string radius of 0.3 - 0.4 fm. The probability 
for a cc pair to pass the central region thus is rather small. This explains the strong 
dependence of the J/"^ results on the string radius. 

We conclude our study by noting that cc dissociation in the prehadronic phase seems 
to be a dominant and likely charmonium absorption process in relativistic heavy ion 
collisions. The qualitative behavior of the proton-nucleus and nucleus-nucleus data can 
be described within this picture. We note that the sudden jump in the J/\E' suppression 
for Pb + Pb suggested by the data for Et ~ 50 GeV in Fig. ^ cannot be described in 
our dynamical model as in Ref. based on string and hadron dynamics (cf. also ||26| 
for a recent overview), because the predicted suppression within our approach is always 
a rather smooth function of atomic number and centrality of the collision. If confirmed 
experimentally, this might indicate a new phase with other degrees of freedom (partons) 



and reaction channels possibly related to a deconfined QGP formation [2S]. 
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Figure 1: Dynamical evolution of a baryonic string. The fragmentation starts after the 
formation time tf. 
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Figure 2: Graphical representation of the hadrons (left) and the strings (right) during 
the high density phase in a central Pb + Pb collision at 160 AGeV in the center of 
mass system. Three time steps are shown at tcm = (4.4, 5.0, 5.6) fm/c (see next figure); 
the axis labels are given in fm. 
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Figure 3: Number of strings and average string length in a central Pb-Pb collision at 
IGOGeV as a function of time in the nucleus-nucleus cms. 
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Figure 4: The J/\l/ survival probability S"^^^ for minimum bias p + A reactions at 
200 GeV assuming a 6 mb cross section for the cc absorption on baryons including cc 
dissociation on strings for two different string radii, 0.1 fm and 0.4 fm in comparison 
to the experimental data from P|. 
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Figure 5: The J/\l/ survival probability S"^/* for S + U at 200 AGeV (upper part) and 
Pb + Pb at 160 AGeV (lower part) as a function of the transverse energy in comparison 
to the experimental data from [Q. The circles in the lower part are the new NA50 data 
from 0. Results are shown for four string radii from Rs = 0.1 fm to Rs = 0.4 fm. 
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